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O ' Abstract 

We address the problem of two pairs of fermions living on an 
arbitrary number of single particle levels of a potential well (mean 
field) and interacting through a pairing force in the framework of 
the Richardson's equations. The associated solutions are classified 
in terms of a number vi, which reduces to the seniority v in the limit 
of a large pairing strength G and yields the number of pairs not de- 
£3 ■ veloping a collective behaviour, their energy remaining finite in the 

G —* oo limit. We express analytically, through the moments of the 
single particle levels distribution, the collective mode energy and the 
two critical values and G~ of the coupling which can exist on a 
single particle level with no pair degeneracy. Notably and G cr , 
when the number of single particle levels goes to infinity, merge into 
the critical coupling of a one pair system G cr (when it exists), which 
is not envisioned by the Richardson theory. In correspondence of 
G cr the system undergoes a transition from a mean field to a pairing 
dominated regime. We finally explore the behaviour of the excitation 
energies, wave functions and pair transfer amplitudes versus G finding 
out that the former, for G > G~, come close to the BCS predictions, 
whereas the latter display a divergence at G cr , signaling the onset of 
a long range off-diagonal order in the system. 

PACS numbers: 21.60.-n, 21.30 Fe, 24.10.Cn 
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1 Introduction 



In a previous paper [1 we examined in detail the problem of the fermionic 
pairing Hamiltonian H pa i r in the simple situation of one pair of nucleons 
coupled to an angular momentum J = living in a set of L single particle 
levels (s.p.L), e.g. in a major shell of the nuclear mean field. This problem 
was solved long ago for the case of one s.p.L hosting n pairs of fermions: yet 
it still presents aspects deserving further investigation in the non-degenerate 
case. 

Actually this is the situation one faces in applying H vair to real systems 
like nuclei and metals, which in fact represent two extreme situations of the 
non-degenerate case: in the former a major shell is typically split into five 
or six s.p.l. of different angular momenta, in the latter the number of non- 
degenerate levels entering into a band corresponds to a significant fraction 
of the Avogadro number. Moreover in a heavy nucleus the number of pairs 
living in a level may be as large as, say, eight while in a metal is one. 

Recently a renewed and widespread interest for H pa i r in the non-degenerate 
frame has flourished in connection with the issues of the physics of ultrasmall 
metallic grains, possibly superconducting [2 , and of Bose-Einstein conden- 
sation [3]. 

When only one pair in L s.p.l. is considered then the eigenstates of the 
pairing hamiltonian fall into two classes: one collective and L — 1 trapped 
(in between the s.p.l.) states. In this connection it was found in jTj that: 

1. the energy E co u of the collective mode is related to the statistical fea- 
tures of the s.p.l. distribution, in the sense that only few moments of 
the latter, beyond of course the strength of the pairing interaction, are 
sufficient to accurately predict E co u\ 

2. the eigenvalues E" r {y — 2, • • • L) of the states trapped in between the 
s.p.l. and hence constrained in the range e v -\ < E\ r < e„, being the 
Ci (i = 1,---L) the single particle energies, belong to two different 
regimes: in the weak coupling regime, for v both small and large, the 
E^ almost coincide with e u (the more so, the smaller the degeneracy 
of the unperturbed s.p.l. is), whereas in the strong coupling regime, 
provided the degeneracy of the s.p.l. increases with u, the E^ r , for small 
u, almost coincide with e v -\. Notably the transition between the two 
situations occurs more and more sharply, namely in correspondence of 
a precise value of the coupling constant G, as L increases. 
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In the present paper we extend the previous analysis considering two 
pairs of fermions living in many levels. First we shall see that in this case the 
eigenstates of the pairing Hamiltonian are conveniently classified in terms of 
a number vi, which provides a measure of the degree of collectivity of the 
states. This number is directly linked to the number Nq, first introduced by 
Gaudin in ref. [3], through the relation V\ = 2Nq and might be considered 
as a sort of "like-seniority" (hence the notation vi) 1 since it reduces to the 
standard seniority v for large G: as pointed out in [5], it has the significance 
of the number of pair energies (see below) which remain finite as G goes to 
infinity. 

Futhermore for the n = 2 system we address the issues of: 

1. expressing the energy E co u of the collective v i = state in terms of 
the statistical features of the s.p.l. distribution, i.e. in terms of the 
moments of the latter; 

2. exploring the existence and providing an analytic expression of the 
critical values of G (we already know from jB] that at most two of 
them might exist on a s.p.l.) which signal the transition of the system 
between two different regimes; 

3. relating these values of G CT to the one previously found in the n = 1 

case; 

4. studying the transition from the weak coupling (where the mean field 
dominates) to the strong coupling regime (where the pairing interac- 
tion dominates) by following the behaviour with G of the pair transfer 
matrix elements between one and two pairs states; 

5. relating the exact and BCS solution in the quite extreme situation of 
few pairs (in fact two) and few s.p.l. (in fact three). 

Anticipating our results, we will find that indeed, under suitable condi- 
tions, two values of G cr exist. They turn out to coincide in the large L limit 
and, notably, they coincide as well with the one found in ref. for the n = 1 
case. Their analytic expression can also be related, as for E co u, to the mo- 
ments of the distribution of the s.p.l. (in fact to the inverse of the latter) and 
they again mark the transition between the two above mentioned regimes. 

1 Note however that, whereas the seniority v counts the fermions coupled to J ^ 0, the 
number vi refers to J = pairs. 
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Indeed for G > G cr the exact solutions of the Richardson equations appear 
to come very close to the Bogoliubov quasi-particle solutions of the reduced 
BCS and, concerning the one pair transfer matrix element, we display their 
divergence in the proximity of G CT . 

Finally we prove that the existence of G cr relates to the trace of H pair . 



2 Classification of the states 

As well-known the eigenvalues and the eigenvectors of the pairing Hamilto- 
nian 

L j v L 

Hpair = Y^ tv ^3vmu^3um v ~ G ^ A^A U , (1) 

u=l m u =—j l/ fj,,u=l 

A» = Ylra^i^ 1 )^'™^^-™^!^ being the J = pair distruction oper- 
ator, are found by solving the Richardson equations [7]. These, for n — 2, 
reduce to the following system of two equations in the unknown E\ and E 2 
(referred to as pair energies): 

1 fi u 2G 



GV^V + tt^V = (2a) 
fa2e^-E 1 E 2 -E l K ' 

fx — 1 



1-GY +— = (2b) 

fa 2e M -E 2 E x - E 2 



or, equivalently, by adding and subtracting the above, 



2-GY—^ GY V = (3a) 

fa2e lx -E 1 fa2e^-E 2 y 1 



-Y ^ +Y n " = o. (3b) 

fa2e fM -E 1 fa2e ll -E 2 E 2 -E 1 y ' 

In (|2I3|) the degeneracy of the single particle energy (s.p.e.) is f2 M and 
G is the strength of the pairing force. Since the Richardson equations deal 
with pairs of fermions (nucleons) coupled to an angular momentum J = 0, 
their eigenvalues, given by E = E\ + E 2 , are those of the zero-seniority 
states (v = 0). Importantly, these eigenvalues display different degrees of 
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collectivity: hence they are conveniently classified in terms of the latter. For 
this purpose, as mentioned in the Introduction, we introduce a number vi, 
that counts in a given state the number of particles prevented to take part 
into the collectivity, not because they are blind to the pairing interaction 
(indeed they are coupled to J = 0), but because they remain trapped in 
between the s.p.L, even in the strong coupling regime. Specifically we shall 
ascribe the value vi= to the fully collective state, Vi—2 to a state set up 
with a trapped pair energy while the other displays a collective behaviour 
and vi—A to the state with two trapped pair energies. 

This classification is equivalent to the one of ref. |3J, where the n pair 
energies are split (for any n) into two classes: those remaining finite as 
G — > oo and the others. 

We shall explore this pattern of states in both the weak and the strong 
coupling regimes, commencing with the former, which is the simpler to deal 
with. 



3 The weak coupling domain 

In the weak coupling limit of course no collective mode develops, hence v\ 
has no significance. Adopting a perturbative treatment we write the pair 
energies E 1 , E 2 as 

Ei = 2e^ + Gx t (4) 
(Gxi being a perturbation). Hence 

2e M - Et Gx % 4" " e M.) " G ^ 
and, expanding in G, the system ^ becomes 

l + ^L- G J2—-^ \ + ^~^Z V + °( G2 ) = ( 6a ) 

l+ %2 '2( e M- e ^) + 2(e w -e w ) + ^ °' ^ 



where the indices \i\ and /x 2 select one configuration out of the unperturbed 
ones. At the lowest order in G, if /i x ^ // 2 , one has = — and the pair 
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energies Ei = 2e w — GVt^. are real. Thus the energy eigenvalue E = E\ + E 2 
becomes 

S = 2(e w +e Ma )-G?(n w +n Pa ) . (7) 
Instead, if /xi = /i 2 , entailing fi Ml > 1, we have to solve the system 

1 + + Q(G 2 ) =0 (8a) 

l + ^-G V — + + 0(G 2 ) =0. (8b) 

x 2 ^ 2 ( e M _ e M 2 ) ^1 - X 2 

A generic solution reads 

Xl,2 = "(ty* - 1) ± V^Mi ~1> ( 9 ) 

showing that i?i and i? 2 are always complex conjugate. Hence the energy 
becomes 

E = 4e Ml - 2G(fi m - 1) , (10) 

which is of course real. 

In comparing (J7J) and ([TU}) with the lowest energy of one pair, which, in 
the weak coupling regime, is [T] 

£ = 2e w -GO w , (11) 

one sees that, while (j2J) corresponds to the sum of two contributions like (JTTj) 
(the two pairs ignore each other), in (fTU|) the Pauli blocking effect appears. 

4 The strong coupling domain 

Here is where V\ has a significance and we deal with the states with v 1 — 
and 2 (the vi = 4 states are of minor physical interest). 

4.1 vi = 

In this Subsection we study the eigenvalue of v\ = in the strong coupling 
limit. Clearly such a state arises from an unperturbed configuration with 
the two pairs in the lowest s.p.l. (if Qi > 1) or in the two lowest s.p.l. We 
already know that in the degenerate case the collective eigenvalue is 

E = 4e - 2G(tt - 1) , (12) 
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e being the energy and Q the degeneracy of the level. The above yields 
the leading order contribution to the energy in the non-degenerate case and 
represents a good estimate when the spreading of the s.p.e. levels is small 
with respect to GQ, as it is natural to expect. 

To show this, following [1 , we introduce the new variables 

Ei-2e 

Xl = (13) 

where now 

^ L L 

e = — ^fi^e^, and fi = ^f^. (14) 
fi=i fi=i 

We further define the variance 



a 



Ij^n^-e)" (15) 



and the skewness 



7 a 3 n 



^-^^(e^-e) 3 (16) 



of the level distribution; moreover we introduce the expansion parameter 



a = W < 17 > 



Then the system (jHJ becomes 



1 / a 2 a 3 \ 1 / a 2 a 3 \ 

2 + — 1 + — + — 7 + ••• +— 1 + — + — 7 + • • • = 

(18a) 



1 A « 2 « 3 \ 1 / a 2 a 3 . , , 
— + — + — T+--- U + — + — T + --- +H = °- 



(18b) 

At the leading order the system is easily solved, yielding 
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Higher order terms can be obtained through a recursive linearisation proce- 
dure. One gets 

(0-1) 9 (0-2) „ (0-4) 

t,, = — ^ L — a- + ct'y- - 

1,2 (0-1) + 7 (0-l) 



vO^T 1 9 O „ O /n/ , . 

± *T T! 2 Q ' (?7^ ±la V^ + 0(a) ' ( 2 °) 

Note that the two pair energies are always complex conjugate. So the 
system's total energy reads 

E-4e (0-1) 2 (0-2) 3 (0-4) 4 

-go" = - 2 ^T~ " 2a (M + 2a 7 (M + } • (21) 

Since the collective energy of one pair of nucleons living in L levels in the 
strong coupling regime is [T] 

^ ! = -l-« 2 + 7« 3 + 0(« 4 ), (22) 

one sees that, when (] > 1, (j2Tj) becomes just twice the value (|2*2*|) and, 
moreover, the imaginary part of x\ and %2 goes to zero as l/\/0. Thus, in 
this limit, the Pauli interaction between the two pairs vanishes, as expected: 
the two pairs behave as two free quasi-bosons condensed in a level whose 
energy is given by (J22|) . 

In Table Q we compare the result (J21|) with the exact one assuming that 
the two pairs live in the first L levels of a harmonic oscillator well with 
frequency uoq. The vi = collective state arises from an unperturbed config- 
uration with a pair in the lowest and a pair in the next to the lowest s.p.l. 
(in units of the oscillator frequency the energy of such a configuration is 8). 
We see that the difference between the two results never exceeds ~ 15 %, 
even for G = G/ftw as low as 0.1. Note also that the energy of the v \ = 
state scales with the size of the well (i.e. it does not depend upon u ). 



4.2 v t = 2 

In the absence of the coupling term (Pauli principle) in the Richardson equa- 
tions (0) the eigenvalues of the v i = 2 states could be simply obtained by 
adding the collective energy Ei carried by one pair and the trapped energy 
E 2 carried by the other pair. 
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G 


a 






EW 




0.1 


0.89 


11.2 


8.21 


6.01 


7.25 


0.2 


0.44 


7.4 


5.91 


5.36 


5.30 


n i 

u.o 


n 9Q 


o.u 


9 fil 




9 ^1 


0.4 


0.22 


-0.2 


-0.95 


-1.08 


-1.11 


0.5 


0.18 


-4. 


-4.60 


-4.69 


-4.70 


0.6 


0.15 


-7.8 


-8.30 


-8.36 


-8.37 


0.7 


0.13 


-11.06 


-12.03 


-12.07 


-12.08 


0.8 


0.11 


-15.4 


-15.77 


-15.81 


-15.81 


0.9 


0.10 


-19.2 


-19.53 


-19.56 


-19.56 


1.0 


0.09 


-23.0 


-23.30 


-23.32 


-23.32 


1.5 


0.06 


-42. 


-42.20 


-42.21 


-42.21 



Table 1: Energies, in units of huoo, of the state Vi — for different values of 
G = G/Hujq at the order 0, 2 and 3 in the expansion parameter a compared 
with the exact ones. The L = 4 s.p.e. levels and the associated degenera- 
cies are those of a 3-dimensional harmonic oscillator, uj being the harmonic 
oscillator constant. 



This situation is recovered in the very strong coupling limit, where all the 
s.p.e. become essentially equal to e and both e and E2 are very small with 
respect to E\. Indeed the first equation of the system (j2J then becomes 

1 + Q 2 - (23) 
G E\ E\ 

yielding 

E 1 = -G(n - 2) , (24) 

namely the energy of the state with two pairs and v = 2 in the one level 
problem. This result, setting a correspondence between states with vi — 2 
and v — 2, connects seniority and "like-seniority" (or the physics of a 'broken' 
and a 'trapped' pair). 

In the non degenerate case, denoting with e[^ and e[ u ^ the first and the 
u-th eigenvalues {y ^ 1) of the one pair equation 
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the coupling (Pauli) term in (J2J) can be approximated for large G as follows 



2G 2G 2G 2 

(26) 



E 2 -E 1 2e u -2e + GQ Q' 

Hence the equations of the system (J2j) decouple and can accordingly be recast 



as 



J__y ^ = (27a) 

("b) 



where 



1 



2 

1 + ^ 



2 

1 - 1 2 



G ( 2 ) G e {u) -E {1) G 



(28a) 



(28b) 



We thus see that in the strong coupling regime the Pauli principle just 
re-scales the coupling constant, differently, however, for the collective and 
the trapped states: indeed the pairing interaction is quenched for the former 
and enhanced for the latter by the Pauli blocking. 

We compare in Table 121 the exact results with the approximate ones 
(eq. 1)27)1) in the h.o. case for L = 3, in the strong coupling regime (G = 5). 
In this example five states exist: of these only two have V\ = 2, namely those 
associated with the unperturbed configurations having either one pair in the 
lowest s.p.l. and one in the highest or both pairs in the second s.p.l. In the 
Table we see that the estimated trapped energy E 2 almost coincide with the 
exact ones, while the results for the collective energy Ei are satisfactory (the 
error being ~ 7%). 

In concluding this Section, we observe that our results in the strong cou- 
pling domain agree with the findings of reference |Hj , where the problem of 
the pairing Hamiltonian for small superconducting grains is studied. Indeed, 
although in jH] the single particle levels are non-degenerate, the conclusions 
of this paper are not affected by this assumption. Hence the results of ref. jS] 
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^ ex act 


jjj exact 


Ei 


E 2 


-33.325 
-33.988 


3.307 
5.720 


-35.711 
-35.711 


3.309 
5.721 



Table 2: Exact and approximate (eq. (|2*7jl ) energies of the states with vi = 2 
for G = 5. The s.p.e. and degeneracies are those of a 3-dimensional harmonic 
oscillator. All energies are in units of hw . The number of levels considered 
is L = 3. 

for two pairs can be recovered from ours by setting Q u = 1 in the formulae 
of this Section. Likewise, our results can be directly derived from ref. [B] 
by setting n=2 and by grouping the single particle levels into L degenerate 
multiplets. 

5 The critical value of G 

In this Section we first search for the critical values of the coupling constant 
G and for their expression in terms of the moments of the s.p.l. distribution 
for the n = 2 system ( Subsection 15. 1|) . Next, in Subsection 15.21 we discuss 
the physics of G cr . Finally, in Subsection 15. 3[ we compare the Richardson 
solution with the BCS one. 

5.1 The analytic expression of G cv 

In the weak coupling regime, when a state evolves from an unperturbed one 
with the two pairs living in the same level, then the pair energies E\ and 
Ei are always complex conjugate. On the other hand when the state evolves 
from an unperturbed one having the two pairs living in two different s.p.l., 
then Ei and E 2 are real. 

By contrast, in the strong coupling regime the pair energies Ei and Ei of 
the v i — state are always complex conjugate. It is then clear that, if the 
degeneracy Qi of the lowest s.p.l. is greater than one, then the pair energies 
Ei and E 2 of the vi = state are complex conjugate in both the weak and 
strong coupling regime and their behaviour with G is smooth. 

On the other hand if Qi — 1, since in the weak coupling limit the two pairs 
must live on different levels, E\ and E 2 are necessarily real in a neighbourhood 
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of the origin, but become complex in the strong coupling regime. Thus a 
singularity in their behaviour as a function of G is bound to occur. 

In this second case it appears natural to surmise that the singularity takes 
place when Ex and E 2 coincide. In fact in this case the Pauli term of the 
Richardson equations diverges, and it must be compensated by a divergence 
in the sum entering into the system (j2J): this can only happen if Ex or E 2 
coincides with an unperturbed eigenvalue. 

To find out the analytic expression of the critical value of G we start from 
a generic solution with two pair energies E\ and E 2 , which evolve with G till 
they coincide at an unperturbed energy 2e„ of a s.p.l. with degeneracy Q u . 
These pair energies must fulfill the G- independent eq. (jUb), that we cast in 
the form 

{E 2 -E l fY- 4t r + 4 = (29) 

1 ' ^ 1 (2e li -E 1 )(2e,-E 2 ) 1 } 

and which allows to express E 2 as a function of E\. To render this relationship 
explicit we set 

Ex = 2e v + x (30a) 
E 2 = 2e v + x<p„{x), (30b) 

where <p u (x) is easily found to read, for vanishing x, 

n u -2± 2y/i - n u 



a real quantity only when Q u = 1. This occurrence is crucial since, in order to 
compensate the divergence arising from the sum appearing in (J2t0, namely 
-GVL v /{2e v - Eh.) = GVtJx, the divergent Pauli term 2G/(E 2 - Ex) = 
2G/[x((p u (x) — 1)] must be real. 

Assuming then Q u = 1, the link between the pair energies is easily found 
by inserting (J30|) into (J29|) and expanding in x. To fourth order in x one 
obtains 

E 2 = 2e u - x ± 2V {2)u x 2 - 4:Vf 2)u x 3 

+ ( ±9Vf 2)u + 2V% )V ± + 0(a* (31) 
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where E\ is hidden in x and the quantities 



are the inverse moments of the level distribution. 

To compute G cr we recast one of the equations of the system (J2J as follows 

i + f-G t 2^-h v ^ = ' (33) 

which is valid in the x — > limit and transparently displays the cancellation 
of the divergences. Clearly the two solutions of are 



Q{v)± = 1 



E gg ± 
2e w -2e„ \ 4^ (2e„ - 2e 
=1(^1/) \ m=i(m/") v m 



(34) 

which actually correspond to two critical values for G. We thus recover the 
results found long ago by Richardson [H] through a somewhat different route. 

The situation is portrayed in fig. ^ where the behaviour of the pair en- 
ergies with G is displayed for a harmonic oscillator well assuming L — 3. 
Clearly in this case the only s.p.l. with Q — 1 is the lowest one, hence the 
pair energies E\ and E 2 , real in the weak coupling limit, coincide at the crit- 
ical point Crcr (their common value being 2ei) and then become complex 
conjugate; the energy of the associated state evolves in the v\ — collective 
mode. By contrast, for the v i = 2 state, arising from the configuration with 
two pairs living on the second level (which is allowed for the harmonic os- 
cillator well), the two pair energies E\ and E 2 are complex conjugate in the 
weak coupling limit, coalesce into the energy 2ei at the critical point G^~ 
and then become real. One of the two solutions remains trapped above 2e%, 
while the other evolves into a collective state: the sum of the two yields the 
energy of the vi = 2 state. 

It is thus plain that for the occurrence of a critical value of G a s.p.l. 
with Q = 1 must exist. Such s.p.l. is the lowest one in a harmonic oscillator, 
hence for this potential well both the ground and the first excited state of a 
n = 2 system carry one critical value of the coupling constant, namely G^ + 
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o 

0.2 0.4 0.6 0.8 1 1.2 

G 

Figure 1: Behaviour with G of the pair energies of the ground {yi = 0) and 
first excited (v t = 2) states obtained as solutions of the Richardson equations 
for a harmonic oscillator well for L = 3: dashed lines denote the real part 
of the solutions (of course coincident) when the solutions are complex and 
solid lines the separate real parts. The pair energies are in units of ftwo. As 
explained in the text that the two critical points refer to the two different 
states. 
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and Gci , respectively. In correspondence of these G the pair energies take 
on the value E\ = E 2 = 2e\. Thus for a n = 2 system two (at most) critical 
points exist on a O = 1 s.p.l. 

Finally, owing to the relevance of the = 1 degeneracy, we consider the 
model of L, for simplicity equally spaced, s.p.l. all having = 1, a situation 
occurring in metals and in deformed nuclei. In this instance two positive G cr 
always exist in the lowest s.p.l. when L > 3 (in fact GiV — > oo for L = 2). 
Moreover a positive implies complex E\ and E% for G < and 

since for small G both the pair energies are real, they should evolve from 
an unperturbed configuration connected with the next higher lying s.p.l., as 
illustrated in fig 121 Numerically we have found, for this model, that two G cr 
appear on the second level when L > 9 and on the third level when L > 16. 
Thus in the O = 1 model for two pairs to form, so to speak, a quartet it is 
necessary that they live on adjacent s.p.l. in the unperturbed configuration. 
Furthermore the more excited the configuration is the more not only G, but L 
as well, should be larger for the merging to occur, a fact clearly reflecting the 
competition between the mean field and the pairing force. Finally observe 
that here, at variance with the finding of ref. G^ +1)+ > this 

is simply because we measure the s.p.e. from the bottom of the well rather 
than from the Fermi surface. 



5.2 The significance of G CI 

In 15.11 we have found that, when a s.p.l. with a degeneracy = 1 exists, 
at least two states of the n = 2 system carry a critical value G CT . It is then 
natural to ask what is the relationship of the latter with the G cr found in 

where a system with n — 1 was addressed. Actually in L G cr marked 
the transition from a mean field dominated regime to the one ruled by the 
pairing interaction (as far as the trapped states are concerned). Moreover, 
in [T] we have proved the existence of G CT for any potential well with s.p.l. 
of increasing degeneracy. 

To connect the values of G CT for the systems with n — 1 and n — 2, 
respectively, we consider again the harmonic oscillator potential. For this 
well, in the n = 1 case, by explicitly computing the variance and the total 
degeneracy of the levels, the critical value of G, in units of huj , is 

Gcr * 4 (35) 



15 



10 




12 3 

G 

Figure 2: Behaviour of the pair energies of three excited states in the case 
of fL = 1. The number of levels considered is L = 9. 
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when L is large. For the n = 2 system first notice that, at large L, in G& 
the moment V\ dominates over P 2 : hence, in this condition, 

G«+ ~ Gg)" . (36) 

Moreover, again for L large enough, 

Pi ^ y (37) 

from where the equality 

G CT {n = 1) = G c > = 2) (38) 

follows in the asymptotic L limit. From this outcome it follows that also 
when n = 2 the relevant dynamical element for G < G& is the mean field, 
whereas for G > G^ is the pairing force, as far as the system's ground 
state is concerned. 

One might ask how this interpretation can be reconciled with the findings 
related to the potential well of equally spaced s.p.l. with Q = 1, for which 
no G cv occurs in the n = 1 system, whereas a Git does exist in the n = 2 
system. 

To understand this occurrence consider the structure of the trace of H pa ir- 
This, in fact, for the n — 1 system and as far as the trapped modes are 
concerned is found, in the strong coupling limit, to read 

S (7) = ^y z ^_ (39) 
L ^ 7 + 2 v ; 

k=2 ' 

where = (E^ — 2efc_i)/(2efc — 2e^i) actually measures the shift of each 
trapped pair energy Ek (k = 2, • • ■ L) from the corresponding unperturbed 
energy 2eu- The above holds for any well with L (large) equally spaced s.p.l. 
with degeneracy 

oc F . (40) 

We thus see that (|3*9*|l yields 



for the harmonic oscillator (7 = 2) and 



E (n*=i) = i (42) 
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for the Q = 1 model (7 = 0). 

Hence for the latter the pairing force is so strongly active among the 
trapped pair energies to prevent the occurrence of a transition from the 
pairing to the mean field dominated physics, whereas the opposite happens 
for the harmonic oscillator well. 

In addition the trace of H pa i r entails that in the Q — 1 model the collective 
mode is weaker than in the harmonic oscillator case: indeed to set it up not 
only a strong G, but a large degeneracy as well, is needed. 

Concerning the n = 2 case in the Q = 1 model, the pairing dominance 
regime, defined by G > G~, is indeed postponed at values of G larger than 
in the harmonic oscillator case. In fact one deduces from (J34j) . in the large 
L limit, the expression 



which, en passant, fixes the domain of validity of the perturbative expansion 
in G PI] for the Q = 1 model. 

One might accordingly conjecture that no transition between different 
regimes will occur in this case among the trapped energies: we are currently 
performing the analysis of the behavior of the Zk versus G for the v 1 = 2 
states to ascertain whether this statement holds true. 

5.3 Comparing the Bogoliubov quasi-particle and the 
Richardson's exact solution 

In this Section we compare the exact solutions of the Richardson equations 
with the BCS solution in terms of Bogoliubov quasi-particles. To this scope 
we self-consistently solve the well-known BCS equations for N fermions living 
in L levels, namely 




y/(e v - Xy + 

n u 2 




(44) 




(45) 




(46) 
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being A = G ^ v u v v u the gap, A the chemical potential and Gvi 
In the BCS framework the excitation energy of a system with seniority v 
is given by the energy of v quasi-particles, each carrying an energy 

E? P = ~ A) 2 + . (47) 

In fig. El we display and compare the exact excitation energies E exc = 
E(vi) — E(g.s.) for a L = 3 harmonic oscillator well and the corresponding 
Bogoliubov's quasi-particles predictions for a v = 2 and a v — 4 state, 
whose excitation energies are 2E® P and 2(£ 1 QP + E® P ), respectively It 




0.5 1 1.5 

G 



Figure 3: Behaviour with G of the excitation energy of a two-pairs system 
in a harmonic oscillator well with L = 3. The exact solutions with vi = 2 
(circles) and v i = 4 (triangles) are compared with the BCS results for v = 2 
(solid line) and v = 4 (dashed line). 

appears that for G larger than the highest critical value (G~ ~ 0.6, see 
fig. 1) both the BCS and the exact excitation energies become linear functions 
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of G and, remarkably, are very close to each other, in particular for Vi = 
2. This result on the one hand shows that the BCS theory is a valid (at 
least for the lowest excitations) approximation of the exact physics even 
when the number of pairs and levels is indeed low, thus strengthening the 
correspondence between seniority and "like-seniority" . On the other hand it 
confirms that the Richardson exact solutions behave like (J47|) in the strong 
G limit, as found by Gaudin 0], who proved that for large L the excitation 
energies are given indeed by sums of terms like (|47|) . It is remarkable that 
this appears to be approximately true already for L = 3. 

Finally we like to remark that the exact energies displayed in fig. |3] cor- 
respond, microscopically, to two particles-two holes (the v\ = 2) and to four 
particles-four holes (the vi = 4) excitations without the breaking of any pair. 
In other words they are associated with the promotion of one or two pairs to 
higher lying s.p.l. 

It is also worth mentioning that while the s.p.e. e u are almost constant in 
G, the quasi-particle energies E® p grow linearly with the latter (for G > G cr ) 
and their spreading in energy is lower than the one among the e u . 

We finally note that it would be interesting to compare the Richardson 
exact result with the energy of the excited state known as "pairing vibra- 
tional" state calculated in the quasi-particle representation. This we leave 
for future work. 

6 The wave functions 

In this Section we examine the wave functions of the states so far discussed. 

As well-known, the v = eigenfunctions of H pa ir for a n pairs system are 
expressed in terms of the collective Grassmann variables $^ (referred to as 
s-quasibosons), according to [T2| 

n 

VU**](m) = ns2M ( 48 ) 

k=l 

where the set of indices (mi, . . . , m n ) labeling the unperturbed configuration 
from where the state develops are collectively denoted by m. Since 

L 

^h-E^h*; (49) 
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can be viewed as the wave function of a pair, (|48|) actually corresponds to 
the wavefunction of n independent pairs. The j3 coefficients are related to 
the eigenvalues Ej. according to: 



^sfSW (50) 



the Ck(m) being normalisation factors. When no confusion arises the index 
(m) will be dropped. 

Since (<I> M , $*) = VL^5^ U , it is convenient to replace the <3>'s and the /3's 
with 

^ = -8= and = v^ep , (si) 

respectively. The normalisation of the state | (m) > then requires 

C k {m) = 1 = . (52) 

In the one pair case the above with k = 1 indeed entails < (m)\(m) >= 1, 
whereas for a two pairs state (|H2*|) can still be used (of course with k = 1 for 
the first and k = 2 for the second pair energy), but the associated state is no 
longer normalised to 1. 

In P] we have investigated the wave functions of a single pair living in 
many levels. We now study the behaviour with G of the wave functions of 
two pairs, our interest being focused on the states undergoing a transition, 
in the sense discussed in Sec. 03 

For sake of illustration we take four s.p.l. of a harmonic oscillator well 
and consider the states v i = (namely the configuration (1711,1712) = (1,2)) 
and one of the v\ — 2 (associated with the configuration (1711,1712) = (2,2)). 
We display m fig. 0] and fig. El respectively, for different values of G, their 
coefficients [3^ as functions of the index fi. 

In the v 1 = case, for G < the coefficients 2) are real, while 

for G > Gcr^ + , they become complex conjugate (hence we display both their 
real and the imaginary part); the opposite occurs in the v 1 = 2 case. 

As it appears in fig. 01 the v 1 = state initially has essentially one pair in 
the first and one in the second s.p.l. As G approaches Gcr = 0.171, both 
pairs almost sit on the first level (thus displaying an apparent violation of 
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the Pauli principle, which actually becomes a true violation at the critical 
point), where however the system cannot live. Finally for G 3> Ga , the 
weight of all the components of /3^(1, 2) and j3^(l, 2) are almost the same, 
thus displaying a collective behaviour. 

The (2,2) vi = 2 state (see fig. |3J) instead starts up with two pairs on 
the second s.p.l. (see the case G — 0.1 in fig. [5]). As G approaches = 
0.287, the two pairs seem to live on the lowest single particle level (again 
forbidden by the Pauli principle). Finally for G ^> G& , all the components 
of /? (1) (2, 2) have almost the same weight (like the components of the collective 
state of one pair) while only the first component of (3^(2,2) is significant 
(like for the trapped state of one pair). 

In conclusion we provide analytic expressions for the coefficients for 
G — » 0, G — » oo and G = G CT . 

The weak coupling limit is immediate and one obtains (7711,7712) = 
5 mifl and /3y (mi,m 2 ) = S m2fl . 

The strong coupling limit is easily handled only when v\ — (indepen- 
dently of fli), when the s.p.e. are small with respect to Ex, E 2 . In fact one 
then gets 

^) = ^(ViPl±i) ■ (53) 

The structure of the wave function at the critical points (assuming the 
transition to occur on the fi — 1 level) is more delicate. To the order x 2 one 
finds 




(57) 



22 



In the above x is connected to G by the relation 
x 2 = T 2 (G - g g ) P (2)1 ± V 'T )2 +0((G- G CT ) 2 ) , (58) 

d/ (2)1 + 4/ (2)1' (3)1 + ' (4)1 

and in ()57|) the double sign refers to the two critical points Gcr' )± . Notice 
that the start linearly in x, i.e., in y/G — G cr , thus displaying a branch 
point in the control parameter G at its critical value. 

7 The transition amplitudes 

In this Section we study the pair transfer amplitudes from the vacuum to a 
one pair state and from a one pair to a two pairs state as a function of G, 
to explore how the transition between the two different regimes previously 
discussed is reflected in the matrix element. 

We thus study the transition amplitude induced by the operator 



3n 



X = E A > E E (-^^IrnA^ , (59) 



which enters into H pa i r , and examine the matrix elements 

< 1 pair | ^ 1 > and < 2 pairs 1 pair > . (60) 

In the fully degenerate case the general expression for these transition 
amplitudes between states with any number n of pairs and seniority v, can 
be obtained by inserting into the pairing spectrum a complete set of states 
\n'v' >. Indeed one has 



< n + l,v\A$\n' ,v' >< n',v'\A\n + l,v > 

n',v' 

= |< n+ l,v\A*\n, v > = ( n + 1_ ^)( 
thus getting for the transition amplitude 



v 
2 



(61) 



< n + l,v\A^\n,v >= (n + 1 - |) (tt-n-^J 



(62) 
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Figure 4: In the L = 4 levels harmonic oscillator well, we display for different 
values of G the coefficients (3^ (or their real and imaginary part) as functions 
of the index /i labeling the s.p.l. for the v i = state. The index k labels the 
pairs. 
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Figure 5: The same as in fig. but for the v\ 
pairs living on the second s.p.e. level. 



2 state, made up of two 
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which, while more transparent, coincides with the one obtaiend in ref. [T3] . 
Note the G-independence of (|62|). which, moreover, is diagonal in the senior- 
ity quantum number. 

For large G we could expect the matrix elements (joTIj) to display an asymp- 
totic behaviour coinciding with (jo^jl . namely 

<l,^ = 0|i t |0> ► Vtt (63) 



^ 2 ^ ^ ^ I V2(n — 1) for^ = 

< ,V " 1 ,Vl> g^oo \y/n^2 for vi = 2 



(64) 



if vi is conserved (or if the above matrix elements are diagonal in 

However, the number Vi is not conserved at finite G: only in the strong 
coupling regime the vanishing of the matrix element between states of differ- 
ent Gaudin numbers occurs. 

Hence we now compute the transition amplitudes at finite G. For this 
purpose we first consider the transition from the vacuum to the one pair 
state, namely 

< M|i f |o >= \/^AM > (65) 

which, using (jHOjl and the Richardson's equations for one pair, namely 



^ _ 1 



^ 2e„ — E G ' 

can be recast as follows 



(66) 



< M|it| >= Ci(m) ^ 9p °» = , (67) 



2e„ - E(m) G 



where the label (m) identifies the unperturbed configuration of the state. 

The strong coupling behaviour of (JBTj) is different in the Vi = and the 
vi = 2 case. For the latter the energy E(m), for G — > oo, remains trapped 
in between two s.p.l.: hence the G-dependent normalisation constant Gi(m) 
will tend to a finite value C£°(m). Accordingly the matrix element (|67|) will 
vanish with G. 

On the contrary, in the v i = case, where the unperturbed pair lives on 
the lowest level, E(l) ~ — Gf2. Hence, for large G, since Gf D (l) ~ Gy/Ti, one 
finds 

< (l)|i f |0 >~ Vti , (68) 
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in accord with (|63|) . Thus at large G, for this matrix element, vi and v 
coalesce. 

Next we study the transition from one to two pairs, namely the matrix 
element of v4J between £>*(m) and *Bi(mi, m2)£>2( m i' m 2)- 
A lengthy calculation yields 



< (mi, m2)\A^\(m) >- 



(2), 



f3^{m)^/n~ v + j3 v {m)JU 



(69) 



the normalisation constant reading 



M 2 = l + 



^[(f(m 1 , m2 f'(m 1)m2 )) , $)(m 1)m2 f»( mi) rn 2 ) 



2 



(70) 



It will be computed in Appendix A. 

Using the Richardson equations to get rid of the sums, we recast (jHHJ) as 
follows 



< (mi,m2)\A^\(m) >- 



1 



2Ci(mi, m^)C 2 {m 1 ,rn 2 )Ci{rn) 



N G[(E l (m ll m2) - E(m)][(E 2 (m 1 ,m 2 ) - E(m)} ' 

(71) 

which is real since E\ and E 2 are either real or complex conjugate. 

We display in fig. |H1 the amplitudes (fTT|) for the transitions from a n = 1 
to a n = 2 system (divided, for obvious convenience, by v^). We consider 
the n = 1 system either in the ground collective state (m = 1) or in the first 
excited trapped state (m = 2). On the other hand the n = 2 system is either 
in the v i = (namely (mi, 1712) = (1,2)) or in the V\ = 2 ((mi,m,2) = (2,2)) 
state. The calculation is performed for a harmonic oscillator well with L = 
4, 10, 20 and 40. 

To understand the behaviour of the curves first consider the weak coupling 
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2 r 0.25 r 




02468 10 02468 10 

GL 2 GL 2 

Figure 6: The transition matrix elements (j71j) divided by for L =4 
(solid), 10 (dashed), 20 (dotted) and 40 (dot-dashed) levels of a h.o. well. 
Upper left panel: transition from m = 1 (collective) state to (m 1 ,m 2 ) = 
(1,2) (vi = 0) state; upper right panel: transition from m = 1 (collective) 
state to (mi,m 2 ) = (2,2) (vi = 2) state; lower left panel: transition from 
m = 2 (trapped) state to (mi,m 2 ) = (1,2) {vi = 0) state; lower right panel: 
transition from m = 2 (trapped) state to (mi, 1712) = (2,2) (vi = 2) state. 
An asterisk denotes the position of the critical points. The straight lines 
represent the asymptotic values for G — > 00. 
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limit, namely (from eq. (|69|) ) 



< (jni,m 2 )\A^\{rn) > 



The above yields: 




(72) 



>/3 


for the transition 


m 


= 1 - 


» (m 1; 


m 2 ) 


= (1,2 





for the transition 


m 


= 1 - 




m 2 ) 


= (2,2 


1 


for the transition 


m 


= 2 - 


» (mi. 


m 2 ) 


= (1,2 


2 


for the transition 


m 


= 2 - 


* (m 1; 


m 2 ) 


= (2,2 



The behaviour around the critical points (which of course depend upon 
L) is quite delicate since both the numerator and the denominator in (|69|) are 
singular when E\ and E 2 tend to 2e\. In fact the numerator in (ffi9~j) vanishes 
like x 2 (see eqs. flrSOj) for the definition of x), but so does the denominator. 
In the end the transition (|69p. for small x, is found to read 



o f _o P ± ^i( m ) P (2) 

<(m 1 ,m 2 )|it| (m)> ^ g 1 ^ ■ (73) 

(t::t /6P ( 4 2) ± 8P (2) Pf 3) + 2P ( 4 4) 

The value of the transition matrix elements at the critical value of G is 
marked by an asterisk in fig. 

As we have seen in Section 15.21 the inverse moments of the level distri- 
bution vanish in the L — > oo limit: hence f)73|) diverges. This occurrence, 
not appearing in fig. El because of the chosen normalisation, relates to the 
ODLRO (off-diagonal long range order) which sets in into a system close to 
a phase transition. 

Specifically the diagonal amplitudes for large L behave according to 



L ^ { V2n 9 (g — G cr ) vi = 

< (mi, m 2 )\A^\{m) > ~ { A _ y (74) 

9 [G - G CI ) vi = 2. 
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The off-diagonal amplitudes behave instead as 

< ( mi ,m 2 )\Am m ) > L ~°° 5 (G - G CT ) (75) 

From the figures it is also clear that the critical value of G increases with 

L and G cr L 2 — > 8 as L — > oo (see eq. f)35jl ). 

Finally we consider the large G limit, ruled by eqs. (f7T|) . (f7H|) and (l8T)|h 
Here first observe that the normalisation constant Af (see eq.()80|) in the 

Appendix) is always finite, reading 

2 - £ for v i = 

(76) 

% fori* = 2 . 

In deducing the above use has been made of the asymptotic expressions for 
Ei and E 2 . 

Considering next the transitions matrix element (JTT)) from a one pair to 
a two pairs state with v\ = or 2, we discuss the two cases: 

i) initial state: n = 1, vi = 0. 

For this, at large G, the energy behaves as — Gil and A/" ~ Gy/Ti. 
Now the n = 2 final state can have: 

a) vi = 0. In this case, using the expressions for the pair energies E\ t 2 

and the normalisation constants Ci^, one finds the asymptotic 
result \/2{VL — 1), which coincides with the degenerate one for 
vanishing v (not vi\); 

b) vi = 2. Here one solution (say Ex) behaves as — GQ, while the other 

(E 2 ) remains finite, as discussed in sec. 14.21 Specifically, one finds 

BU ~ -G (1) Q ~ 

Gcff ^" fi + 2 

and E 2 much lower than E\ and of the energy of the initial state. 
The normalisation constants are d = G { £Vn and G 2 °° (this one 
turns out to be G-independent and ~ 1). Using the above and 
the expression (|75|l for A/", one then gets for the matrix element 
the asymptotic value 

g 2 °° rir~ 

G V 0-2 ' 



30 



which is vanishing as G 1 in the strong coupling limit, as it should, 
since it entails a variation Avi = 2 of "like-seniority" . 



ii) initial state: n = 1, vi = 2. 

For this the asymptotic values of C and of the energy (denoted by C°° 
and E°°) are finite: C°° is of order 1 and E°° is of the order of the 
unperturbed levels. Again, the n = 2 final state may have: 

a) vi = 0. In this case E°° is negligible with respect to E\ and E 2 , 
whose asymptotic behaviour is known. The matrix element is then 
found to read 

V2C°° 



G,yQ{Q-i) ' 

Again a transition Avi = 2 entails a behaviour of order G~ x , but 
in this case the extra factor ~ Q^ 1 induces a faster decrease of the 
transition matrix element, as clearly apparent in fig. (lower left 
panel) ; 

b) vi = 2. Calling E%° (and likewise C%°) the G —>■ oo values of the 
pair energy (and of the normalisation constant) of the trapped 
pair, we find the following expressions for the asymptotic matrix 
element 

2C 2 00 C 00 



y/VT^2G{E ao - Ef) 

Now two different possibilities occur: if E°° and E%° are not both 
confined in the range (2e M _i, 2e M ) then, in the strong coupling 
regime, their difference remains finite and the same occurs for the 
normalisation coefficients C°° and C%°. Hence the matrix element 
is ruled by the factor 1/G and thus vanishes. Otherwise the differ- 
ence E°° — E%° tends to vanish. In fact, let E(G) be the solution 
of the one pair equation Then the trapped pair energy E 2 , 

in the strong coupling limit, solves eq. (|27b ). and, as a conse- 
quence, one has E 2 = E(G)J), with G e g given by eq. (|2*5jl . Thus, 
inserting E 2 = E — y into eq. (|2~7b ) and expanding in y, one finds, 
using ((HSJ), 

2\C°°\ 2 

y = -GDT- 
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Hence the G — > oo limit of the transition matrix element becomes 

where the ratio / C°° can be shown to be given by 

Cf G_ _ _ 2 

^eff 

Therefore the asymptotic matrix element reads y/Q — 2, namely 
the value of the degenerate case (|55)l. 

8 Conclusions 

It is well-known that the pairing Hamiltonian, a reduced version of the BCS 
model of superconductivity, is solved by the Richardson's equations as far as 
the states of zero seniority v are concerned. 

In this paper we have explored the solutions of these equations in the 
simple case of only two pairs of fermions, but living in any number of s.p.l. 
of any potential well. In this framework a further classification of the states in 
terms of Vi, which accounts for their degree of collectivity, appears convenient. 
It is elegant that V[ and v coincide for G — > oo, whereas at finite G one (vi) 
classifies the states according to the same pattern of the other (v) without 
breaking any pair. At small G, of course, no need for vi is felt. 

Concerning the structure of the solutions we have naturally searched, and 
obtained, an expression for the energy of the collective v i — v — mode in 
terms of the statistical features of the s.p.l. distribution, namely in terms of 
the moments of the latter. The same we have achieved for the two critical 
values of the coupling which mark the transition from the physics of the mean 
field to the one of the pairing force. 

According to the general theory of Richardson for a two pairs system at 
most two critical values (G+ and G~) can exist on an unperturbed s.p.l. 
provided the pair degeneracy of the latter is one. Remarkably and G CT 
are found to merge in the L — ► oo limit and here coincide with the G CT , 
previously found in pQ, for a system of just one pair living on any number 
of s.p.l. Note that in the latter G cr should not exist according to 

Richardson, but we have found that in fact it does by following the behaviour 
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of the trapped modes with G. Its existence essentially rests on the trace of 
H P air or, equivalently, on the balance between the collective and the trapped 
modes, a balance reflecting not only the strength G of the pairing force, but, 
critically, the degeneracies of the s.p.l. In the n = 2 case the latter markedly 
affect the critical values of G, which, in turn, signal the onset of the validity 
of the BCS framework. Indeed, for example, we have found that the G CT of 
the Q — 1 model is much larger than the one of the harmonic oscillator well. 
Whether this will entail the absence in the Q = 1 model of a transition in 
the trapped modes for the n = 2 system, as it happens for the n = 1 system, 
is an issue we are currently investigating. 

Be as it may, we have explored in detail the transition associated to 
and G~ T by following the behaviour with G of both the excitation energies of 
the vi = 2 and Vi — 4 states and of the pair transfer matrix elements, finding 
for the latter a most conspicuous enhancement (the more so, the larger L) in 
the proximity of G cr , clearly reflecting the onset of an ODLRO in the system. 
Concerning the excitation energies we have found that they closely approach 
the BCS predictions at large G especially for the case of two quasi-particle 
(vi = 2) excitations, notwithstanding that the latter, in the BCS picture, 
amounts to break a pair, an occurrence never happening in the Richardson 
frame. 

Finally we have found that the BCS ground state of the n = 2 system in 
the simple L = 3 harmonic oscillator model has A/d ~ 2.5 in the proximity 
of the critical values, being d the spacing among the s.p.l. (actually A/d is 
growing linearly with G). The Anderson's criterion for superconductivity 
is thus fulfilled. Indeed in the ground state our system can be viewed as four 
fermions (a "quartet") sitting on the same energy level. 

A Appendix 

In this Appendix we shortly provide analytic expressions for the normali- 
sation constants, whenever this is possible. When both the pair energies 
are real no simple formulas for the normalisation coefficients Ci^tti, n) are 
available, whereas when E{ = E 2 one can use eq. (|2*5j) to get 

(77) 

The above allows to simplify the expressions for the normalisation con- 



Ei — Ei 
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stant distinguishing, however, whether the solutions are real or complex con- 
jugate. 

In the former case a tedious calculations provides 



M= I 2[Cf(m,n) + Cl(m,n)} ^ 
y [E 2 (m, n) - Ei(m, n)] 2 



(78) 



still in terms of the unspecified normalisation constants C\ and C 2 . 

In the latter case the normalisation is known, but, to get A/", we need to 
define 

C k (m,n) = —= 1 , k =1,2 (79) 

^ {2e^-E k {m,n)f 

where, unlike in (JH2|) . no absolute values appear. Thus the C&(m, n) are 
complex and (7| = C\. With the help of (|75jl. and using (f77f). one then gets 



34 



References 

[1] M. B. Barbara, R. Cenni, A. Molinari and M. R. Quaglia, Phys. Rev. C 
66 (2002) 034310. 

[2] G. Sierra, J. Dukelsky, G. G. Dussel, J. von Delft and F. Braun, Phys. 
Rev. B 61 (2000) 11890. 

[3] J. Dukelsky, C. Esebbag and P. Schuck, Phys. Rev. Lett. 87 (2001) 
066403. 

[4] M. Gaudin, "Modeles Exactement Resolus", les Editions de Physique, 
France, 1995 

[5] J. M. Roman, G. Sierra and J. Dukelsky, Phys. Rev. B 67 (2003) 064510. 

[6] R. W. Richardson, Jour. Math. Phys. 6 (1965) 1034. 

[7] R. W. Richardson and N. Sherman, Nucl. Phys. 52 (1964) 221. 

[8] E. A. Yuzbashyan, A. A. Baytin and B. L. Altshuler, Phys. Rev. B 68 
(2003) 214509. 

[9] M. Hasegawa and S. Takaki, Phys. Rev. C 35 (1987) 1508. 

[10] J. M. Roman, G. Sierra and J. Dukelsky, Nucl. Phys. B 64 (2002) 483. 

[11] M. Schechter, Y. Imry, Y. Levinson and J. von Delft, Phys. Rev. B 63 
(2001) 214518. 

[12] M. B. Barbara, A. Molinari, F. Palumbo and M. R. Quaglia, Phys. Lett. 
B 476 (2000) 477. 

[13] D. Bes and R. A. Broglia, Lezioni di Varenna International School of 
Physics "Enrico Fermi", Course LXIX, edited by A. Bohr and R.A. 
Broglia, Varenna, 1977, p. 59. 

[14] P. W. Anderson, J. Phys. Chem. Solids 11 (1959) 28. 



35 



